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1. The Langevin equation and NSPT 

A typical task for Lattice Perturbation Theory (LPT) is the determination of renormalization 
factors or the separation of non-perturbative effects from observables related to confinement. The 
gluon and ghost propagators belong to this class. ^ For that higher-loop calculations are needed. 
Unfortunately, the diagrammatic approach is much more involved than in the continuum case. It is 
unlikely to go beyond the two-loop level in near future. 

A promising alternative to diagrammtic LPT has been proven to be Numerical Stochastic Per- 
turbation Theory (NSPT). For a review see Ref. [Q]. It makes it possible to obtain higher-loop 
results without computing vast numbers of Feynman diagrams. Various applications of NSPT have 
been described in the past, see e.g. [^, 0, ^ ||]. Here we report on first NSPT studies of the Yang- 
Mills gluon and ghost propagator in Landau gauge. For other preparatory results see also [Q, ^J. 

NSPT is based on stochastic quantization, i.e., a lattice variant of the Langevin equation. This 
equation describes the evolution of 4D fields with respect to a fictitious time t under the influence of 
Gaussian random noise. In the limit t —^oo the whole set of gauge links U is distributed according 
to Gibbs measure exp(— Sg[U]). In our simulations we use the Wilson gauge action Sg[U], but an 
improved action would not be much more complicated. 

Other than in usual Langevin simulation, in NSPT the gauge link fields are expanded in powers 
of the bare coupling g oc jS^'/^. Discretizing the Langevin time as t = nt, rescaling T = e/p and 
using the Euler scheme, a set of coupled equations emerges corresponding to different orders in the 
coupling constant. 

From the resulting fields the Green functions of interest can be constructed order by order in 
LPT. In NSPT the algebra-valued gauge potential A^^^ are related to the gauge lattice link fields 

Ux,^i by 

A.,^=Iog[/,^. (1.1) 

Its expansion is given in the form 

^x,^i 2^P ' x,i>- ~ ■ y^--'-) 

I>0 

Each simultaneous Langevin update of the expansion coefficients A^||j is augmented by a stochastic 
gauge-fixing step and by subtracting zero modes from A''' as described in Ref. |@]. 

In this talk we present perturbative contributions to the lattice gluon and ghost propagators 
in minimal Landau gauge. Before taking measurements, this gauge is achieved by iterative gauge 
transformations. We use a perturbatively expanded version of the Fourier-accelerated gauge-fixing 
method 

2. The gluon and ghost propagators in NSPT 

The lattice gluon propagator D"^y (q) is the Fourier transform of the gluon two-point function, 
i.e., the expectation value 

DfM = (A^(fc)At(-fc)) = 5'"^D,y{q) . (2.1) 
' See the talk by A. Maas at this conference nil] trying to identify the non-perturbative content. 
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{k) is the Fourier transform of A" ^, and q denotes the physical discrete momentum correspond- 
ing to the integers G (— L^/2,L^/2] , 



'7m (^P 



- sm 
a \ Ln 



nki. 



(2.2) 



In NSPT the different loop orders n (even orders in (5 '/^) are constructed directly from gauge 
fields A;;^'') (A:) 



2n+l 



1=1 



1^ {k)Ay 



■b,{2n+2-l) 



i-k) 



(2.3) 



,{0) 



Already the tree level term D^y results from quantum fluctuations of gauge fields with / = 1. Terms 
with non-integer n = 1/2,3/2, .. . in (23) have to vanish numerically. Motivated by the structure 
of the continuum propagator in Landau gauge we consider the so called gluon dressing function of 
different loop orders n 



(2.4) 



The ghost propagator is nothing but the inverse of the Faddeev-Popov (FP) operator M which 
can be constructed in Landau gauge using the lattice covariant and left partial derivatives. Since 
the progagator is color diagonal, G"^(^) = d"''G{q), it is obtained as the color trace 
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Gi4) = j;^^{Tr.,iM-\k))^. 



(2.5) 



In (2.5) M ^{k) is the Fourier transform of the inverse FP operator in real space. The perturbative 
expansion is based on the mapping {A^'Ij} {M^'^} {[M^^]^'^}. Since M is expanded in terms 
ofM^ (containing A^'^), a recursive inversion is possible: 



[M 



-il(o) 



-1 



, [M 



-nC) 



/-I 



-11 oo 



(2.6) 



The momentum-space ghost propagator at ?i-loop order is obtained from even orders 1 = 2nofM ' , 



sandwiching (2.6) between the plane-wave vectors \k,c) 



\k,c). 



(2.7) 



Again non-integer n (odd /) orders have to vanish within numerical precision. We present our 
results in terms of the ghost dressing function: 



(2.8) 



Note that the perturbative construction of M in terms of the gauge fields A (and therefore the cor- 
responding ghost propagator) differ from the definition adopted in most Monte Carlo calculations. 
For each chosen momentum (^1,^2,^3)^4) and different colors c of the plane wave \k,c) the prop- 
agator has to be calculated individually. This makes the measurement expensive. 
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3. Some numerical results 

The configuration sequence of all A^'^ created at finite e can be used to measure the pertur- 
batively expanded observables. Already at finite e the non-integer n contributions to the dressing 
functions have to vanish. We have studied the step size limit £ — s- working at £ = 0.07 . . . 0.01. 
In order to make contact with standard infinite-volume LPT at vanishing lattice spacing the limits 
L — > OO and — > have to be performed additionally. 

In Fig. [l] we present different orders of the gluon dressing function at the smallest time step £ = 
0.01 and lattice size 10^. The loop contributions are labelled by integer n, dressing contributions 
with noninteger n vanish. As an example of the ghost propagator we show in Fig. |2| the one- and 
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Figure 1: Measured gluon dressing function Jgi{q) vs. q^. Left: Separate loop contributions J^^\q) vs. q^ 
(for inequivalent lattice momentum 4-tuples) at L = 10 and £ = 0.01; right: the vanishing coefficient of the 
contribution j3^"/^. 

'^(1) '^(2) "(1/2=3/2) 

two-loop results Jl and Jl for the dressing function together with the vanishing Jl 
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Figure 2: Measured ghost dressing function /gh [q) vs. q^ (for inequivalent lattice momentum 4-tuples close 
to the diagonal) for lattice sizes L = 6, . . . , 16 and for the time step e = 0.01 . Left: The one-loop and two-loop 
contributions; right: the vanishing coefficient of the contribution ^ fi^^l^. 



For each set of inequivalent lattice momenta (^1,^:2,^3,^4) we have performed the extrapo- 
lation to zero time step. Performing that limit, the individual loop contributions to the dressing 
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functions are available for all lattice momenta at a given lattice size. They can be compared with 
Monte Carlo results on finite lattices. 

The perturbative dressing functions summed to loop order n^nax are calculated for a given 
lattice coupling p as follows: 



4l/gh(^'"max) — o„''gl/gh 
n=0/lP 



(3.1) 



In Figs. ^ we present the summed dressing functions at fixed jS as function of the number Wmax of 
loops. Remarkably, all loop contributions are of the same sign, such that we expect that these 
summed-up dressing functions represent a sequence of lower bounds for the total perturbative 
function for all momenta q^. The variation with the lattice coupling is shown in Figs. ^ The 
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Figure 3: Summed gluon (left) and ghost (right) dressing functions for j3 = 6.0 up to four (left) and three 
(right) loops for a set of momentum 4-tuples. 
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Figure 4: Summed perturbative gluon (left) and ghost (right) dressing functions at different /3 . 



large loop contributions to the gluon propagator that we find can be interpreted as resulting from 
lattice tadpole contributions. 

4. Summary 



In NSPT gauge link variables and gauge potentials are naturally related via ( 1 . 1 ) different from 
most nonperturbative numerical implementations. Using this definition, we have performed higher 
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loop calculations of the lattice gluon and ghost propagators in Landau gauge to make predictions for 
the perturbative content of the propagators as function of the lattice momenta taking the hypercubic 
group into account. 

To compare with Monte Carlo data and in this way to find out the nonperturbative contributions 
of the propagators, the logarithmic definition of the gauge fields and the corresponding FP operator 
have to be implemented. Work in that direction is in progress in the Humboldt University group. 

A first attempt to estimate the unknown two-loop contribution to the lattice gluon and ghost 
propagator in Landau gauge in the limits of infinite volume and aq —^0 has been presented recently 
on other occasions and can be found in the Proceedings p^ . 
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